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Метод решения гамильтоновой задачи 
коммивояжера 


Предлагается двухэтапный метод поиска решения гамильтоновой задачи коммивояжера, который 
либо находит решение поставленной задачи, либо корректно устанавливает, что задача неразрешима. 
Разработанный метод имеет значительно меньшую потребность в вычислительных ресурсах, чем 
известные алгоритмы. 


Введение 


В статье формулируются признаки негамильтоновости графа, опираясь на 
которые еще на первом этапе решения задачи можно определить, разрешима ли она. 
Это дает возможность на втором этапе искать оптимальное решение. 

Сформулируем гамильтонову задачу коммивояжера (ГЗК). 

Задан связный граф Н = (Г, И) без петель и кратных рёбер, Г — множество 


вершин, 


/|=и, И -— множество рёбер. Ребру {1 ЛИ приписан вес (стоимость) 


а; Е 7, Ы=Ьп, 25 — множество неотрицательных чисел. Симметричная матрица 


+ ле : ь а 
весов [4], где Ч; Е 2, если {1 ЛЕЦ, и 4; =® иначе, {#], 4; =, 1=1,п, 
полностью определяет взвешенный граф Н = (Г, 0). 

Граф Н = (Г, И) гамильтонов, если он содержит гамильтонов цикл, то есть 


простой цикл, проходящий по всем вершинам Г точно один раз. Гамильтонову циклу 2 
поставлена в соответствие стоимость (2), равная сумме весов его рёбер. 


ГЗК состоит в нахождении гамильтонова цикла 2* минимальной стоимости 0(2*). 

Поставленная задача №Р — трудна [1]. Известные методы её решения представ- 
лены схемами организации полного перебора всех циклов в графе Н [2-5]. Практи- 
ческая реализация этих методов проблематична даже с применением самых мощных 
компьютерных систем. 

ГЗК характеризуется тем, что она не всегда разрешима. Поэтому единственным 
методом поиска гамильтонова цикла минимальной стоимости 2* является метод, 
который либо строит 2*, если ГЗК разрешима, либо корректно устанавливает, что 
граф Н = (У, И) негамильтонов. 


Основной результат 


Поиск решения будем выполнять в два этапа. На первом этапе проверяются 
достаточные условия негамильтоновости графа Н = (У, И). ГЗК не разрешима, если 
выполнено хотя бы одно из достаточных условий негамильтоновости графа. 

Очевидно, что если граф Н содержит концевые вершины, то ГЗК не имеет 
решений. Менее очевидным является условие негамильтоновости графа, содержащего 
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точки сочленения. Точкой сочленения называется вершина, при удалении которой 
вместе с инцидентными ей рёбрами, образуется несвязный граф [6], [7]. Иначе говоря, 
если у точка сочленения связного графа Н, то граф Н -— ъ несвязен (рис. 1). 


ть 


а) 6) 
Рисунок | —а) Граф Н с точками сочленения В и]: 
6) компоненты связности, как результат удаления вершин В и / 


Утверждение. Если граф Н = (ТУ, И) содержит точки сочленения, то он 
негамильтонов. 

Доказательство. Вершина у связного графа является точкой сочленения тогда 
и только тогда, когда существуют две другие вершины х и у, такие, что любая цепь 


между х и у содержит вершину у. Но в гамильтоновом цикле, кроме цепи (х, ..., у, ..., у) 
должна быть цепь (у, ..., х), которая не проходит через вершину у. 

Следствие. Пусть связный граф Н = (У, И) содержит цикл 2 = (р, 4, ..., Ьр), где 
6,>2, 0, =...=0, =2. Тогда граф Н негамильтонов. 


Доказательство. Рассмотрим пару вершин Г и /], ГЕФ, ] #2. Так как граф 
связен, то любая цепь, соединяющая вершины / и /, проходит через вершину р. 


Следовательно, вершина р является точкой сочленения. 
На рис. 2 представлен граф Н, удовлетворяющий условиям следствия. 


Рисунок 2 — Граф Н содержит цикл (3, 2, 1, 4, 3), где д. =4, 


9) =д, =9, =2 и поэтому негамильтонов 


Распознавание точек сочленения в связном графе Н = (Г, И) выполняется с 
трудоёмкостью ОИ] [6]. 

На втором этапе ищется решение ГЗК 2* для графа, не содержащего концевых 
вершин и точек сочленения. 

Алгоритм поиска 2* построен в соответствии с классической схемой метода 
ветвей и границ, вызывающей процедуру решения задачи о назначениях (ЗН) для 


вычисления оценок снизу величины Ь(2* 4]. Вместе с тем, он имеет присущие 
У 


только ему особенности, которые в процессе ветвления позволяют определить, какие 
подзадачи ГЗК не имеют решения. 
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Рассмотрим матрицу весов [4 | ГЗК. Чтобы найти оценку снизу искомой ве- 
* п 
личины О (2 — ‚ требуется для этой матрицы решить задачу о назначениях (ЗН). Но ЗН 


с матрицей стоимостей [4 | содержит часть элементов а; = © и может не иметь 
. п 


решения. Следовательно, для вычисления нижних оценок в процессе поиска 2* 
требуется алгоритм, который или корректно находит решение ЗН, или строго 
фиксирует, что ЗН неразрешима. Так работает модификация алгоритма Кана- 
Мункреса, предложенная в [4]. 

Алгоритм Кана-Мункреса решает ЗН на максимум в предположении, что все 
элементы 4, =®, 1=7. Входом модифицированного алгоритма служит матрица 


[4 | ‚ где 4, =а-4,, если а, #0, и а = о иначе, 4 = — максимальный эле- 
я 


мент матрицы Га, . Тогда, если существует решение ЗН на максимум для матрицы 
С п 


[4 | ‚ то оно является решением ЗН на минимум для матрицы [4 | . Стоимости ре- 
п з п 


шения С (л) и С’(л) соответственно для матриц [4, |. и [4 | связаны равенством 
С(л)=па-С'(л). 

Модифицированный алгоритм Кана-Мункреса ищет решение ЗН в двудольном 

графе К(Х, У, Е), |Х|=У|=пв, Е=20 


‚ отвечающем матрице |4, | ‚ в котором вер- 
. п 


шина х, е Х соединена с вершиной у; ЕТУ ребром (х я с весом 4 (х у,)= а, #0. 
ЗН разрешима, если в графе К (Хх ‚У,Е ) построено совершенное паросочетание л с 


максимальной суммой весов рёбер. ЗН неразрешима, если граф К (Хх ии: ) не 


содержит совершенных паросочетаний. 
Детальное описание модификации алгоритма Кана-Мункреса представлено в [4]. 
Её основной конструкцией является известная процедура поиска совершенного паро- 


сочетания в двудольном не взвешенном графе Н (х ,У,Е ) ‚© | = |у Е | = 20 


ненная средствами, устанавливающими при каком количестве элементов а; =—© и 


‚ ДОПОЛ- 


> 


их расположении в матрице [4 


В ЗН неразрешима [4], [6]. Как алгоритм Кана- 


Мункреса, так и его модификация характеризуется трудобмкостью, оцениваемой 
величиной 0(»*) [6]. 


Алгоритм поиска 2* выполняется по схеме метода ветвей и границ, предложенной 
в [2] для решения симметричной ЗК. Эта схема в комбинации с модификацией 
алгоритма Кана-Мункреса применима и для решения ГЗК. 

Пусть построено совершенное паросочетание л. Оно доставляет целевому функ- 
ционалу ЗН минимум С (л) ‚ принимаемый за нижнюю оценку стоимости искомого 
цикла 2*. Рассматривая паросочетание л как перестановку столбцов матрицы 
стоимостей, представим его цикловым разложением, то есть в виде совокупности 
непересекающихся циклов. Цикловое разложение перестановки л и оценка С(л) 


образуют корень дерева ветвлений. Перестановка л, представленная единственным 
циклом, является решением 2* ГЗК. Его стоимость равна С (л) : 
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В общем случае решение ЗН содержит несколько непересекающихся циклов. 
Выберем из них цикл а = (и и: у) с наименьшим числом рёбер. Удалим 


все решения ЗН, которые содержат цикл о’, не потеряв при этом в исходной матрице 
ни одного допустимого решения 2. Поскольку хотя бы одно из А рёбер (у,.›,), 
(У›,уз),... (Ум) не должно входить в 2, то множество всех решений ЗН можно 
представить разбиением на А подмножеств. ЗН с исходной матрицей [ау |. обозна- 
чим Ё. Тогда А, разбивается на К подзадач А, Р...., В,. Каждой из них взаимно 


однозначно соответствует ребро цикла о’, вес которого полагается в матрице [4 | 
п 


равным со, а все остальные веса остаются без изменений. В матрице [а, | ЗН на 
ы п 
максимум этому ребру приписывается вес —о. Тогда, если существует решение 
ГЗК, то оно принадлежит множеству решений какой-либо из подзадач В, Р,,..., 
В. ‚ представленных вершинами дерева ветвлений, исходящими из вершины А. 
В каждой подзадаче Р, 1=1, К, можно добиться исключения не только решений, 
содержащих цикл о’, но и решений, включающих циклы с вершинами из множества 


9 = у р - С этой целью в матрице стоимостей подзадачи Р, полученной из 


В присвоением элементу 4,„„, 


= 1 А, У=и ‚ веса со, положим Ч,» = © для всех 
УЕЗ\ {у} . В соответствующей матрице стоимостей подзадачи ЗН на максимум каж- 
дый элемент Ч,» получает вес —©. 

Для ЗН на максимум, соответствующей подзадаче Р, {= 1А ‚ построим моди- 
фикацией алгоритма Кана-Мункреса перестановку л;, если Р разрешима, или 
определим, что она не имеет решений. Если подзадача Р неразрешима, то соот- 
ветствующая ей вершина в дереве ветвлений не имеет допустимого продолжения. 
Предположим, что из К подзадач Р разрешимы А подзадач Е. ЕК 


К 


5 = 1, К ‚› То есть построены оптимальные перестановки п; и вычислены доставляе- 


мые ими значения С(л., | . Очевидно, стоимость искомого гамильтонова цикла "А 


можно ограничить снизу величиной 
С = тт С (т, | |1 Е{!, 2,..., ®}, 5 Г} 


Рассмотрим подзадачу Р,, для которой С (л.) =С. Если решение л, является 


гамильтоновым циклом, то оно будет и решением ГЗК. В противном случае 
перестановка Пл образует несколько непересекающихся циклов. Тогда вершина В, 


дерева решений объявляется вершиной ветвления [2]. Задача Р, разбивается на под- 


задачи, решения которых не содержат цикла минимальной длины о’ из разложения 
перестановки т. Ч. В решениях подзадач исключены также все циклы, порожденные 


на множестве вершин о. Выполнив модификацию алгоритма Кана-Мункреса для 
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каждой полученной подзадачи, определим А› разрешимых подзадач Р,, Че\1, 2, 


... К}, 5=Ь К. . Теперь текущей нижней границей стоимости оптимального гамиль- 
тонова цикла 2* является величина 


@ = тя [ттт С (ло, | |4 Е {1, ря ®}, 5 1%}, тт С (т, )! 


Ге 1, ыы к} ; 5Е11, Ра К} 5Ь #4}}, 


которой соответствует задача т Если все подзадачи, полученные из задачи ыы 
неразрешимы, то вершине ветвления дерева решений соответствует задача Р, с 


текущей оценкой 
С = тт С(л, ) [2=112,... 52, „вр 9] 


Дальнейшее ветвление для вершины В осуществляется точно так же, как и для 
вершины Р,. Поиск решения ГЗК завершается в одном из двух случаев. В первом 


случае алгоритм находит гамильтонов цикл 2* ГЗК, который является решением ЗН 
с текущим значением С. Во втором случае алгоритм устанавливает, что все 
конечные вершины дерева решений не поддаются дальнейшему ветвлению, и, 
следовательно, ГЗК неразрешима. 


Пример 

Для графа Н = (У, И), изображенного на рис. 3, и входной матрицы стоимостей 
|4; | выполним поиск решения ГЗК. Граф Н не имеет висячих вершин и точек 
сочленения. Поэтому нельзя утверждать, что граф Н негамильтонов. 


Для матрицы [4 |. выполним поиск решения ЗН А, применив 


модифицированный алгоритм Кана-Мункреса [4]. 


1.2 34 615 
о х 9 5 10 6 1 4 
о ® © © [1 4 9 4 
Эх ® 4 хх ® 2 ® 
5 © 4 ® хх П 9 ® 


10 Шо о ох 12 о © 


11 9 2 9 © © © © 


© чм фр - 


4 4 © © © © © © 


Рисунок 3 — Граф Н = (Т, И) и матрица стоимостей его рёбер 


634 «Искусственный интеллект» 32008 


Метод решения гамильтоновой задачи коммивояжера Г 


Алгоритм ищет решение ЗН на максимум для матрицы. 


1 2. 3 4 ть 6 7 8 
1 —со —со 3 7 2 6 1 8 
2 со —оо —со —ю 1 8 Э 8 
о —оо —оо 8 —ю© —ю 10 —о 
|4 р ] за 4 |7 —со 8 —оо —ю 1 3 —со 
5 2 1 со со о 0 —© —со 
6 |6 8 —о 0 со со со 
7 1 3 10 3 со со со со 
8 18 8 со со со со со со 


Оптимальным решением ЗН К как на максимум, так и на минимум является 
перестановка л = (4, 8, 7, 1, 6, 5, 3, 2), С’(л)= 50. Минимум целевого функционала ЗН 
В, С(л)= 46 ограничивает снизу стоимость искомого решения 2* ГЗК. Цикловое раз- 
ложение перестановки л имеет следующий вид: о\ = (1, 4, 1), о, = (2, 8, 2), о.= (3, 7, 3), 
од = (5, 6, 5). Каждый цикл разложения содержит два ребра. Поэтому для ветвления 
выберем любой из четырёх, например, о. = (3, 7, 3). На рис. 4 представлено дерево 
ветвлений, отображающих процесс поиска гамильтонова цикла 2* С(л). Все результаты 


вычислений, формирующие дерево ветвлений, содержатся в табл. 1. 


С(л)=46 


Рисунок 4 — Дерево ветвлений для ГЗК с числовыми данными примера 1 


ЗН А в результате ветвления порождает две задачи А и Р› на максимум с 


' 


< ! 
матрицами стоимостей, полученными из @ | В матрице [а, |. для задачи Р на —0 


заменяется элемент 457 =10, а для задачи Р, — элемент 4..=10. Модифицированный 
алгоритм Кана-Мункреса находит решения л| и л>› этих задач и определяет стоимости 
полученных решений С(л!) = С(л.) = 49. Для ветвления можно выбрать любую из 


подзадач с одинаковыми оценками. Выберем задачу Р. Перестановку л| исчерпывают 
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два цикла, из которых цикл (3, 4, 7, 3) имеет минимальную длину. Исключение 
решений ЗН, содержащих этот цикл и всех циклов с вершинами из множества {3, 4, 7}, 
даёт три задачи А‚, Н,, И, на максимум и соответствующие им матрицы стоимостей. 


В первой матрице 4 =4;. =—9, во второй 4. = 4; = 441 = —0, в третьей 
4; = 413 = 414 =—%. Все три задачи разрешимы. Стоимости их решений равны 
С (м, )=57, С (ль, =: (^, ) = 58. На текущий момент все концевые вершины 
дерева ветвлений Р,, А, Н,, А, активны. Так как С = тт {С(л.), С(т, }, С( т, ) | 


С (*,, ) = тт 449, 57, 52, 58} = 49, вершина Р, оказывается вершиной ветвления. 


Таблица 1 — Результаты вычислений, формирующие дерево ветвлений 


Цикловое разложение решения 
ЗН 


(1,4, 1) (2,8, 2) (3, 7,3) (5, 6, 5) 
(1,5, 6,2, 8, 1) (3, 4, 7, 3) 

(1, 5, 6,2, 8, 1) (3,7, 4, 3) 

(1,4, 7, 3, 1) (0, 8,2) (5, 6, 5) 


(1, 8, 2, 7, 3, 4, 1) (5, 6, 5) 


ЗН | Решение ЗН 

В, | л=(4, 8,7, 1,6, 5, 3,2), С(л)=46 
В |л.=(5, 8, 4,7,6,2,3, 1), С(л1) =49 
Ь |1›=(5,8,7,3,6,2, 4,1), С(л>)= 49 
В |7 =44, 8, 1,7,6,5,3,2) С” )=57 


1 


В, | ль= (8,7, 4,1, 6, 5,3,2), С(л,)=52 
В; | ль=(,7,4,3,6,5,2, 1), С(л,)=58 


(1:80) (2,71,2)(3,4,3)6,6;5) 


Р; п = (8, и. Ио 6, р 4, 2), С ло) 


(1, 8, 2, 7, 4, 3, 1) (5, 6, 5) 


2 


(1, 4, 3, 7, 2, 8, 1) (5, 6, 5) 


(т, ) = 63 
В, | ло,= (4,8, 7,3, 6, 5,2, 1), С(т, )=52 
В, | лз=(3, 8,7, 1,6, 5,4,2), С(л,)=57  |(1,3,7,4, 1) (2, 8,2) (5,6, 5) 
Ви | Льл= (8, 7,4, 6, 1,5,3,2), С(л»)=56 |(1,8,2, 7,3, 4, 6,5, 1) 
В>› | ЗН неразрешима = 
Ро | Поз= (5, 8, 7,3, 6, 4,2, 1), С(л»}=56 |(1,5,6,4,3,7,2,8, 1) 


Р,>› | ЗН неразрешима Е 


Задача Р, разбивается на три подзадачи Р,, Р,,, Р», со стоимостями решений 


С( т, ) = 63, С(л,, ) =52, С( т, ) =57 . Из текущего множества активных вершин {В : 
ВВ В» В; дерева ветвлений выберем вершину с минимальной оценкой 


С =тт 457, 52, 58, 63, 52, 57} = 52. Равноправными претендентами на вершину 
ветвления являются вершины А,, Р,›. Для дальнейшего разбиения выберем подзадачу 


А›. В циклическом разложении перестановки л\, цикл (5, 6, 5) имеет минимальную 
длину. Поэтому задача А, порождает две подзадачи А, и В››. Модифицированный 
алгоритм Кана-Мункреса строит решение ЗН В›| в виде гамильтонова цикла (1, 8, 2, 7, 3, 


4, 6, 5, 1) стоимостью С (711 ) = 56 и устанавливает, что задача В›› не имеет решений. 
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Метод решения гамильтоновой задачи коммивояжера Г 


С построением гамильтонова цикла стоимостью 56 отпадает необходимость 
ветвления всех активных вершин с не меньшими оценками. Остаётся единственная ЗН 
Ь, ‚ стоимость решения которой равна 52. Циклическое разложение перестановки л›, 
включает цикл (5, 6, 5), вызывающий две новые подзадачи В.1 и Р,››. ЗН Р.>› не имеет 
решения, а решением ЗН Р,›, является гамильтонов цикл (1, 5, 6, 4, 3, 7, 2, 8, 1), у 
которого суммарный вес рёбер С(л»>\) равен 56. 


Стоимость построенных гамильтоновых циклов меньше, чем нижняя граница в 
любой концевой вершине дерева ветвлений. Таким образом, оптимальным решением ГЗК 


являются ЦИКЛЫ 27 = Л]! = (1, 8, 2, 7, 3, 4, 6, 5, 1), 22 ЕО = (4956: 3, 72.8, 1). 


Заключение 


Метод реализован программно на языке С#. Вычислительный эксперимент 
проведен на машине Зеегоп 1,8 Гц с оперативной памятью 2 Гб. Генерировалось не 
менее 100 матриц размерности 40, 60, 80, 100, которые заполнялись на 75 % 
значениями из интервала 1 — 50, а на 25 % — ©. Полученные во время эксперимента 
результаты содержатся в табл. 2. 


Таблица 2 — Результаты, полученные во время эксперимента 


Размерность матриц 40 60 80 100 
Минимальное время решения, с 10,25 15,78 23,27 39,18 сек 
Максимальное время решения, мин | 4,58 26,22 65,26 121,75 мин 
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1.В. Гаращенко, А.В. Морозов, А.В. Панйиев 

Метод розв?язання гам!льтоновот задач! комвояжера 

Пропонуеться двоетапний метод пошуку розв’язання гампльтоново! задач! комвояжера, який або 
знаходить розв’язання поставлено! задач, або коректно встановлюе, що задача не мае розв’язання. 
Розроблений метод мае значно меншу потребу в обчислювальних ресурсах, ниж в1дом1 алгоритми. 


Т.И. Сагазйсйепко, А.Г. Могосоу, А.Т. Раш5йеу 

Тье Мево4 Гог Зо!ушо о? Нат Иоп?$ Соттегса! ТгауеПег РгоМет 

ТБе 5\о-$азе те@о4 1$ оеге4 ог зо[уш5 Фе Нат Июп соттегсла] {тауеПег рго ет у\/ысЬ ейфег Нпа$ а 
зомНоп ог сотесу зрес1Иез Фаё фазК ш чпзоуаЫе. ПОеуе]оре@ тео Баз з1ютийсапИу 1еззег 
сотршайопа| ро\ег гедитетеп($ ап Кпо\уп апаюзчез. 
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